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Abstract: 

We obtain the operator product expansion of the self-energy in the O(N) non-linear 
cr-model to all orders in the coupling and the large momentum, and to next-to-leading 
order in 1/N. In the light of this result we discuss recent suggestions that there may 
be additional power corrections from short distances, associated with defining the cou- 
pling constant non-perturbatively. The non-linear a-model provides no evidence for such 
'non-standard' power corrections. We also find that the OPE converges for sufficiently 
large external momentum, presumably because there are no multi-particle thresholds at 
arbitrarily high energies in the 1/N expansion. 



dk 2 F(k,Q)^ (1) 



1. Short-distance observables in QCD, characterized by a large scale Q ^> Aqcd, can 
be predicted perturbatively in a s (Q) < 1. In some cases non-perturbative power cor- 
rections in Aqcd/Q can be incorporated as well, using the operator product expansion 
(OPE) 0, [|. Alternatively, one can deduce the structure of non-perturbative effects 
by investigating the infrared (IR) contributions to loop integrals in perturbation theory. 
This the basis of the so-called renormalon method [H] . Either way, for example for current 
two-point functions % Jdxe iqx (0\T(j(x)j(0))\0), where j is a bilinear quark current, one 
finds that the leading power correction scales as (Aqcd/Q) 4 and is related to the gluon 
condensate. Both methods parametrize power corrections that arise from long-distance 
contributions to the observable. In general one cannot exclude non-perturbative effects 
from short distances. The only known dynamical mechanism for such corrections is due 
to instantons of small size ^> Aqcd H- They scale as (Aqcd/Q) 11 ~ 27V///3 and are 
strongly suppressed at large Q. 

In recent work Grunberg || and Akhoury and Zakharov |J suggested that one should 
expect quite generically power corrections of order (Aq CD /Q) 2 from short distances on 
grounds that the perturbative QCD coupling is usually not specified to an accuracy 
better than this. To make their point let us consider the integral 

ajjc) 
k 2 

which is typical of renormalon calculations of diagrams with one gluon line. When the 
loop momentum k is small, and if F(k, Q) behaves as {k/Q) p for small k, one obtains a 
contribution of order (A/Q) p to the integral. This term stands for the power correction 
that is taken into account by the OPE or the renormalon method. In this context the 
definition of the coupling a s has played no role. It is usually considered to be the one- 
loop perturbative running coupling, which has a Landau pole at k 2 = Aq CD . However, 
physical quantities do not have a Landau pole and a non-perturbative definition of a s (k), 
valid at all k, should not have a Landau pole either. Such a non-perturbative coupling 
can itself be expanded in a power series, when k is large: 

^ = R 1 (1-2IK2 \ +■■■+ COnSt • ^ + • • • • ( 2 ) 

-fa ln(A; 2 /A QCD ) k 2 

The dots denote further perturbative and power correction terms. If the second term 
is present, and if one insists on extending (|l|), derived in the context of perturbation 
theory, to a non-perturbative running coupling, one obtains a power correction of order 
Aq CD /Q 2 to the integral in ([TJ) from the region of large k ~ Q. In [|5], |5| these additional 
power corrections are interpreted as a potential new source of non-perturbative effects, 
not accounted for by the OPE, but also not in contradiction with the OPE, because they 
arise from short distances (large k). The concept of an effective coupling is not unique 
and the existence of a particular 1/k 2 term in (|2]) is not compelling. Independent of 
the precise form of power corrections, however, the argument given in |5], [J suggests the 
existence of unconventional power corrections at some order. 
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The present work is motivated by our attempt to understand this argument, and its 
relation to the discussion of a related issue in J7|, in the context of the two-dimensional 
non- linear O(N) a-model. We calculate the asymptotic expansion for the self-energy of 
the a particle to next-to-leading order in 1/N and to all powers in the large momentum. 
The approximation contains all non-trivial elements of the OPE construction in QCD, 
except for the absence of multi-particle production. In particular, it allows us to address 
the question of non-perturbative corrections to the effective coupling and potential power 
corrections from short distances. Our conclusions can be summarized as follows: 

(a) There are power corrections from the region k ~ Q, but they are included in the 
standard OPE construction and correspond to soft subgraphs of the graphs that define 
the effective coupling non-perturbatively. Despite the fact that k ~ Q, these power 
corrections are therefore of infrared origin. 

(b) The fact that the perturbative coupling has a Landau pole does not give rise to 
additional power corrections. The Landau pole disappears and the correct analyticity 
structure is recovered only after summation of the OPE. The issue of defining an effec- 
tive coupling non-perturbatively is not related to the existence of power corrections to 
particular physical observables. In fact, it turns out that the OPE is naturally expressed 
in terms of the perturbative coupling (the first term on the right-hand side of (^)). 

(c) The OPE for the self-energy is convergent for \Q 2 \ > 9m 2 in the complex Q 2 
plane, where m is the dynamically generated mass of the a particle, analogous to Aqcd 
in QCD. The radius of convergence is determined by the onset of the 3-particle cut. The 
convergence of the OPE is probably a consequence of the 1/N expansion, because to 
the order we are working at most three particles can be produced. More generally, we 
expect that the OPE of the self-energy is convergent for \Q 2 \ > (2k + l) 2 m 2 if the 1/N 
expansion is truncated at the order k. Thus, for finite N, we expect the OPE to be 
divergent for any Q 2 . Because of the restriction on the number of produced particles, 
the 1/N expansion in the cr-model is not suited for studies of the analytic continuation 
of the OPE to Minkowskian momenta and the validity of "parton-hadron duality" , and 
we do not address both questions in this work. 

We emphasize that while the present calculation clarifies the argument of [|5], |5| it 
does not preclude the existence of power corrections from short-distances due to a yet 
unknown dynamical mechanism other than small instantons, see || for a particular sug- 
gestion. Up to this date, however, all (known) short-distance power corrections are of 
semiclassical origin and unrelated to integrals of the type (|l|) and the issue of defining 
the coupling constant. 



The nonlinear O(N) a- model was considered repeatedly as a toy model for the OPE 



10, 11]. In particular, the non-perturbative separation of long-distance and short- 



distance contributions is well understood in this model, both in cut-off schemes |L0 
and in dimensional regularization ||. In the second factorization scheme, the coefficient 
functions have non-Borel-summable perturbative expansions due to infrared (IR) renor- 
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malons. The prescription dependence that follows from defining the series is exactly 
compensated by the prescription dependence of defining the power-divergent matrix el- 
ements of higher-dimension operators. In the subsequent calculation we observe this 
compensation to all orders in the OPE. However, in this Letter we are mainly concerned 
with power corrections not related to IR renormalons, in particular the effects of a cou- 
pling redefinition on the power series expansion, and with analytic properties of the 
power series truncated to finite order. 

The construction of the 1/N expansion and the Feynman rules are detailed in P, |10[ . 
The N a-particles have massQ 

m 2 = /i 2 e -l/ 9 M (3) 

where \l is a renormalization scale. In addition to the standard massive propagator 
for the scalar cr-particles, the 1/N expansion involves the propagator of the auxiliary 
"ct-particle" 



D(k) = -4tt (k + 2m 2 ) g eS (k) = -AnJk 2 (k 2 + 4m 2 ) In 



Vk 2 + Am 2 + Vk 2 ' 



\/k 2 + 4m 2 - Vk 2 



(4) 



In order to establish a connection with (|1|), we interpret g e s(k) defined through the a- 
propagator as above as a non-perturbative effective coupling. The factor k 2 + 2m 2 is 
separated such that g e s{k) can be viewed as the sum of the chain of a-particle bubble 
graphs with momentum-dependent vertices 



= l/ a ,-I( t ,m) = » + * O + 9 ° OO + 



where an UV regularization is implied, go is the bare coupling and 

Tii \ q(q + k) , , 

^ m) = l d q [ q 2 + m 2 ][( q + k) 2 + m 2 y (5) 

The coupling g e s{k) can be expanded at large momenta generating a power series 

2 

g cS (k)=g(k)~2g 2 (k)— + .... (6) 

Following the terminology of 0, g c s(k) naturally splits into a perturbative and non- 
perturbative part, g e s(k) = g{k) + Sg(k), where the perturbative part g{k) is given by 
the first term on the right-hand side alone. 



1 To avoid factors of 4-7T, we rescale the coupling by a factor 47r. With this convention g satisfies 

ar \ 2 d 9 2 
PK9) = I" U-o = ~9 ■ 



Note that this /3-function is exact in leading order of the 1/N expansion. In particular, g[k) has a 
Landau pole at k 2 = to 2 . 
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(a) (b) 

Figure 1: cr-self-energy diagrams at order 1/N. 

It is important to note that the power corrections Sg(k) to the running coupling are 
of IR origin. Diagrammatically, the reason for this is that even if the momentum that 
flows into the chain of bubbles is large, one of the lines in any one or more of the bubbles 
can still be soft. 

D(k) does not have a Landau pole, but a cut at —k 2 > Am 2 . Because k 2 + 2m 2 
is factored out, g e s{k) has a 'kinematic' singularity at — k 2 = 2m 2 in addition. Note, 
however, that whenever the power expansion of the effective coupling is truncated, there 
is an additional Landau singularity at k 2 = m 2 . 

3. Consider the propagator of a er-particle with momentum p 2 ^> m 2 . To leading order 
in 1/N it is equal to (p 2 + m 2 ) -1 and we can interpret its expansion at large momenta 
as the operator product expansion 

I oo / 2 \ "i oo 

(^A'^E^ ^£c<»V)<o,.><°>, (7) 

V n=0 \ V J y n=0 

with C^\p 2 ) = —l/p 2n and (O n )^ = m 2n . The OPE interpretation implies that the 
terms identified as operator matrix elements originate from large distances. To see this, 
one has to start from the perturbative phase of the a model and construct the OPE 
explicitly; details can be found in |TUfl . We shall take the OPE interpretation of the 
large-iV result for granted, and concentrate on the first subleading order in the large 
iV-expansion, formulated in terms of massive a-fields already. 

To 1/N accuracy the cr-propagator involves the self-energy correction given by the 
two diagrams shown in Fig. [I], where wavy lines denote the a-propagator (|J). The self- 
energy is quadratically ultraviolet divergent. We define the renormalized self-energy by 
zero-momentum subtractions 

= E(p) - E(0) -P 2 ^m\ p2=0 (8) 

and drop the superscript "ren" in what follows. The tadpole diagram in Fig. |I|b is then 
subtracted completely, and the self-energy is given by 

If k 2 + 2i7i 2 

= ^iV J Sk (p + ky + 9cS ^ + subtractions - ( 9 ) 
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We have to verify that the asymptotic expansion of the propagator can be written in the 
factorized form 



{p 2 + m 2 )- 1 - (p 2 + m 2 )- 2 £(» = 

1 CO 

= -2 E [c£V) + (V^V)] [(O n )W + (l/A0(On> (1) ] + 0(l/iV 2 ) (10) 



n=0 



where C^ 1 - ) (p 2 ) and (O n )^ correspond to 1/N corrections to coefficient functions and 
matrix elements, respectively. Note that the explicit construction of the operators is 
not necessary and is complicated in the case at hand by the existence of infinitely many 
operators of the same dimension. The OPE statement means that the contributions from 
small and large distances to (§) can be factorized. In the present case this also means 
that all contributions from large distances are associated with operator matrix elements. 

Returning to the integral @, long-distance contributions originate from the IR sensi- 
tive regions k 2 ~ m 2 and (p+k) 2 ~ m 2 . The OPE is obtained by expanding the integrand 
of ([)]) in small quantities given a particular loop momentum region. The following three 
loop momentum regions have to be considered separately: 

(i) k 3> m, p + k ^> m: Both lines are far off-shell and both the cr-propagator and the 
effective coupling can be expanded in m 2 . The diagram effectively contracts to a 
point and the result can be interpreted as a 1/N correction to the short-distance 
coefficients times operator matrix elements evaluated at leading order in 1/N. 

(ii) k ~ m, p + k ^> m: The cr-propagator can be expanded, but the effective coupling 
cannot. Thus the a propagator is contracted, but the a-propagator is not; the 
resulting diagram can be viewed as a 1/iV-correction to the vacuum expectation 
value of the product of two a-fields (0|a(x)a(0)|0) with small separation x ~ 1/p. 
Expansion in x yields as the leading contribution a 1/iV-correction to the matrix 
element (0|« 2 |0) times a leading-order short-distance coefficient. (This matrix 
element was evaluated explicitly in |], IHJ.) Note that from power counting it 
follows that this region contributes to (|9]) first at order p 2 ■ m 4 /p 4 , i.e. with a 
l/p 4 -suppression similar to the gluon condensate in QCD. (The argument does not 
hold for the contribution of the subtraction terms.) 

(iii) k m, p + k ~ m: The effective coupling can be expanded, but the cr-propagator 
cannot be expanded. This case looks symmetric to the previous one but the in- 
terpretation is different: contracting the a-propagator one gets a leading contribu- 
tion at large iV to the vacuum expectation value for the product of two cr-fields 
(0|cr a (x)cr a (0)|0) = l/g{x' 1 ) ■ [1 + m 2 x 2 /A + . . .] (no sum over 'a'). Note the in- 
verse coupling as a consequence of the normalization of the "length" of the cr-field. 
The result can therefore be interpreted as a 1/iV-correction to the short-distance 
coefficient times a leading large N operator matrix elements of cr fields. Note that 
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a contribution of this type is usually not considered in the renormalon method. 
However, it could be accounted for as well. In the QCD analogue of our problem, 
this would correspond to renormalon contributions to the coefficient function of a 
quark condensate. 

We emphasize that all three contributions are part of the usual OPE construction and 
that (i,iii) include the "short-distance" power corrections discussed in As already 

mentioned, these power corrections are in fact also of long-distance origin. 

For the actual calculation of the expansion of @ one could follow the above procedure 
and calculate the contribution from each region separately. It is more convenient not to 
perform this split-up and to calculate the expansion directly. To this end we write the 
inverse logarithm in D(k) as 

■A-V* 



dt 



where A = (1 + Am 2 /k 



2\l/2 



A + l 

Subsequently, we use the Mellin-Barnes representation 



(11) 



\A-V 




IA + 1. 





A 

where the kernel K(s, t) is given by 

r(t + s)T(-l -2s) 



ds 
2^1 



K{s,t) 




(12) 



K(s,t) 



r(t - i - s) 



1 2{t + s) (t + s)(t + s + l) 



t-8-1 (t-S-l)(t-s) 



(13) 



After these manipulations all integrations can be done and we arrive at the following 
"Borel representation" for the expansion in powers (and logarithms) of m 2 /p 2 : 



-t/g(p) 



1 



ng[t\V (14) 



Explicit expressions for the functions H^[t] are given in the appendix. 

All dependence on m 2 /p 2 is either explicit or implicit in g(p) = l/ln(p 2 /m 2 ). 
We first explain the origin of the various terms in ([Rj). 

The functions F^ 1 [t] and [t] correspond to the contributions of regions (i,iii) dis- 
cussed above. Because they are multiplied by the factor e - */ 9 ^, they can be interpreted 
as Borel transforms of the perturbative expansions of short-distance coefficients times 
matrix elements that do not depend on g(p). The separation of regions (i) and (iii) is 
not unique. In region (i) the expansion of the cr-propagator in m 2 /(p + k) 2 leads to IR 
divergences, which require an intermediate regularization. The result from region (i) 
then takes the form 

I 1 /'ii ii.i 1 j-i (15) 
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where /i is the intermediate factorization scale. On the other hand, the expansion of 
g e ff in region (iii) leads to increasingly UV divergent integrals, for which the same 
intermediate regularization is required. The result from region (iii) then takes the form 

FW[t] ln-^ + G«[*]. (16) 

Both contributions combine to the square bracket in (|14|). The intermediate (and ar- 
bitrary) scale fi drops out and the logarithms combine to the inverse running coupling. 
Recall that region (iii) is a long-distance contribution to the integral @, but it is multi- 
plied by a perturbative expansion, because large momentum flows through the effective 
coupling. After combining (i) with (iii) as in (0), both F^[t}/g(p) and G^[t] have 
long- and short- distance contributions. 

The function H$[t] originates from region (ii) discussed above and therefore should 
be identified with the 1/N correction to operator matrix elements. It cannot be writ- 
ten as a series expansion in g(jp) and has only long-distance contributions. No extra 
regularization is needed to separate (ii) from (i,iii). 

Nevertheless, the separation of (ii) from (i,iii) is not unique. The functions F^[t], 
[t] , [t] have singularities (typically simple poles) in the complex t plane at integer 
values t = ±fc, k = 1,2,.... These are just the usual ultraviolet and IR renormalon 
singularities. (The poles at t — correspond to the usual logarithmic UV divergences 
and are cancelled by renormalization.) Using the expressions collected in the appendix, 
one can verify that all IR renormalon singularities at positive values of t cancel. The 
cancellation of a particular singularity at t = to occurs between G^[t] and H^ +t ^[t] 
and thus involves a cancellation between a short-distance coefficient and an operator 
matrix element over different orders in the power expansion. This is an explicit all-order 
verification of the cancellation mechanism discussed in ||. As a consequence of the 
singularities in individual terms of the sum over n, the summation and the integration 
over t cannot be interchanged, unless the integration contour is shifted slightly above (or 
below) the real axis. Only after such a definition can one truncate the OPE. Obviously 
the prescription dependence cancels in the end. 

As an explicit example, consider the leading asymptotic behaviour of £(p) as p — > oo, 
related to the expansion of F^ n =^[t] and G^=^[t] around t -> and to #<£ =0) [t]. We 
find 

,2 



(17) 



\ng(p) + 1.887537- 2 g(jp) + Y,^g n+1 (p) 

n=l 

with factorially divergent coefficients 

<r n = n!{[l + (-l)»]C(n + l)-2}. (18) 

Because of the Riemann zeta-function this gives rise to an infinite series of renormalon 
poles. 
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Returning to fll4"|), we note that the OPE is naturally expressed in terms of the per- 
turbative coupling g(p). The reason is that the asymptotic expansion is mathematically 
an expansion in powers and logarithms of m 2 /p 2 , and absorbing logarithms only into the 
coupling leaves all powers of m 2 /p 2 explicit. This implies that at every finite order in the 
OPE there is a Landau singularity at p 2 = m 2 , which disappears, when the expansion is 
summed. 

One can avoid unphysical Landau singularities order by order in the asymptotic 
expansion, if one eliminates g(p) in favour of g e ff(p), which amounts to an all-order re- 
arrangement of the OPE. However, it is important to stress that there is no advantage 
to doing this: although the Landau singularity never appears, the predictive power is 
not enlarged. One would still have to sum the entire OPE to be able to predict the 
correlation function accurately in the region p 2 ~ m 2 . In general, (but not in the case of 
g e ff(p)), an inadequate choice of the non-perturbatively defined coupling may obscure 
the cancellation of renormalons poles, or even introduce power corrections into the OPE 
which are not natural for the physical quantity in question, but which appear in order 
to cancel power corrections implicit in the coupling definition (see the discussion in 0). 

In the present example all power corrections do indeed originate from long distances 
and can be accounted for in the OPE framework. The calculation clarifies that the fact 
that the perturbative coupling has a Landau pole, or that a non-perturbative coupling 
has itself power corrections does not provide sufficient motivation for introducing ad hoc 
power corrections beyond the OPE. Moreover, all known sources of short-distance power 
corrections are related to semiclassical effects. Semi-classical effects are connected with 
multi-particle emission, a phenomenon which is certainly beyond the approximation in 
which one assumes the emission of one quantum gluon through an effective coupling. It 
therefore seems that to address the question of whether the OPE is complete, one would 
have to abandon the effective coupling representation. 

A short comment is in order on ultraviolet renormalons, which lead to factorially 
divergent, sign-alternating series in the cr-model (as in QCD). This divergence limits the 
accuracy of a purely perturbative approach and the best possible approximation has an 
error that scales as a power in m 2 /p 2 . It is sometimes argued (see for example ||) that 
UV renormalons require adding power corrections (of short-distance nature) to the OPE. 
In our example, the series expansion of the leading term (|1|) as well as all others lead 
to UV renormalon singularities at negative t. Being located at negative t, these singu- 
larities do not lead to ambiguities in the individual terms of (|T4]) . Since the OPE (|TJ]) 
reproduces the exact result (as one can check numerically), we also conclude that UV 
renormalons do not necessitate introducing additional power corrections, as expected on 
general grounds for sign-alternating series. 

4. Given the OPE to all orders, it is interesting to ask about the convergence/divergence 
properties of the OPE itself (rather than the perturbative expansion that multiplies each 
power correction). For example, in |T2|] a toy heavy-quark correlation function is used 
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to demonstrate that in this case the heavy quark expansion is factorially divergent. 

For the self-energy in the cx-model it is straightforward to derive from the formulae 
collected in the appendix that the contribution of large n S> t (for fixed t) to the sum in 
(|l4l) is equal to 



95/4 m 4 




irp 2 



- \n(9m 2 /p 2 ) + ip(t) - ^(1 - 1) - r(t)r(i - t) 



(19) 

It follows that the sum converges pointwise in t for \p 2 \ > 9m 2 and uniformly on any 
interval [0, to]- The contribution to from t > to can be made arbitrarily small by 
increasing to as can be seen from inserting ( |ilf ) with lower integration limit to into (f|). 
Hence the domain of convergence of the OPE (i.e. the series after t-integration) is given 
by \p 2 \ > 9m 2 . 

The radius of convergence seems to be related to the fact that the self-energy graph 
has a discontinuity for — p 2 > 9m 2 , related to the fact that any cut contains three cx- 
particles of mass m. We suggest that the convergence of the OPE is a consequence of 
the 1/N expansion. At every finite order of the 1/N expansion, a cut can contain only 
a finite number of particles. We expect that the OPE of the self-energy is convergent 
for \p 2 \ > (2k + l) 2 m 2 if the 1/N expansion is truncated at the order k. For finite N, 
when there is no restriction on the number of particles in a cut, we expect that the 
OPE diverges for any value of p 2 in agreement with the conclusion of fll2 |. It may be 
interesting to pursue further the connection between the divergence of the OPE and 
multi-particle production. 

Another consequence of the absense of multi-particle thresholds is that, viewed as a 
function of g, the self-energy has no other singularities in g than cuts stretching from 
g — to g — l/(ln9 + inn) (n a non-zero integer) in the complex coupling plane. This 
differs from QCD and the cx-model at finite N, where the analyticity domain has zero 
opening angle fl3[], because in the presence of multi-particle thresholds up to arbitrarily 
high energies the singular points accumulate at g = on curves with zero opening angle 
at g = 0. 

5. In conclusion, the expansion of the self-energy in the cx-model provides the first 
analytic example of an OPE to all orders in m 2 /p 2 in the non-trivial situation that each 
term in the expansion is multiplied by an infinite series in the coupling that contains UV 
and IR renormalons. We discussed the cancellation of IR renormalon ambiguities in this 
expansion. We found that the OPE is convergent in the 1/A-expansion and converges 
to the exact result without the need for additional power corrections. 

In [TJ[] an unexpected A 2 /p 2 - correct ion (where A is the QCD scale) is reported in 
the OPE of the plaquette expectation value. The present analysis does not offer an 
explanation of this fact. One difference is that the lattice calculation is done at finite 
UV cut-off h.jjy = a -1 with p and Kjjv being identified. One then obtains additional 
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power corrections from higher- dimension operators to the action. We can mimic this 
effect by adding the dimension-4 operator C/Afj V ■ a 2 d 2 a to the cx-model action. The 
new vertex contributes an additional factor k 2 /Afj V to the integral @. For kjjy ~ p this 
results in an order 1 contribution to regions (i,iii) and an order m 6 /p 6 contribution to 
region (ii). The analogous argument for Yang-Mills theory shows that one cannot obtain 
a A 2 /p 2 -contribution in this way. The only loophole is that for Aj/y ~ p one should 
consider the infinite series of higher-dimension operators. If this series is divergent, this 
may induce further power correction terms in the lattice theory at finite cut-off not 
present in the continuum theory. 

Acknowledgements. We thank G. Grunberg and G. Marchesini for useful discussions. 
The work by N. K. was partially supported by NORDITA through the Baltic Fellowship 
program funded by the Nordic Council of Ministers. 

Appendix. We collect explicit expressions for the coefficients in the asymptotic expan- 
sion of the self-energy in (|14|). For n = 0: 

= 1, 

G p 0) M = | + r ^ T -V'(l + t)-V'(2-t)- 2 7 £, 
1 



H$jj[t) = - + B 1 (t). 



(20) 



For n — 1: 



1 
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G®[t] = - 7 + 2t 2 -4t + l + (l-t 2 )^(l + t)+^(2-t) + 2 7i? 

1 



- + — + — — - B (t) 
t t-l t+1 W 



(21) 



where ip(x) is the digamma-function and B and Bi are scheme-dependent subtraction 
constants which are analytic in t: 



B (t) = 3 J(t, 1), B x {t) = 6J(t, 2) - 6 J(t, 3) + 6J(1 + t, 3) - 7J(1 + t, 2) 



J(t, m) 



3^2 



1, m, t + 1 



t + 1 

For generic n > 2 we obtain: 



1 _i_ * — — |— — 

~ 2 ' 2 ~ 2 



fc=0 



k\k\ 



(22) 
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n-2 



H^[t] = Y,K np (n-k-2,t) 

n 

G^[t] = J2K p (n-k,t) 



(n-k-l)t 



k\k\ 

n-k-l+t)l 



(23) 



fc=0 



k\k\ 



il>(t) - -t) + 2ip(k + 1) - 2^(n + t - 1) 



where (z)k = T(z + k)/Y{z) and 
1 T(2t + 2n- 1) 



K p (n,t) 



n\ Y{2t + n- 1) 



+ 



n(n - 1) 



#n P (M) = 



T(-t, 1 +t,2n + 3) 
r(t + n + 1,3 + n-t) 



2t + n-l (2t + n)(2t + ra - 1) 



( 2 n + 2-t (t-2-n)(t-l-n) 
t + n + 1 (t + n + l)(t + n + 2)_ 



(24) 



In the last formula we used a shorthand notation T(a, b, c, . . .) for the product of T- 
functions with the respective arguments. 
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